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Abstract 

Preliminary group classification for a class of generalized inviscid Burger's equa- 
tions in the general form u t + g(x, u)u x = f(x, u) is given and additional equivalence 
transformations are found. Adduced results complete and essentially generalize re- 
cent works on the subject . A number of new interesting nonlinear invariant models 
which have non-trivial invariance algebras are obtained. The result of the work is 
a wide class of equations summarized in table form. 
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1 Introduction 

The theory of Lie groups has greatly influenced diverse branches of mathematics and 
physics. The main tool of the theory, Sophus Lies infinitesimal method [5], establishes 
connection between continuous transformation groups and algebras of their infinitesimal 
generators. The method leads to many techniques of great significance in studying the 
group-invariant solutions and conservation laws of differential equations [21 [HI [10]. The 
approach to the classification of partial differential equations which we propound is, in 
fact, a synthesis of Lies infinitesimal method, the use of equivalence transformations and 
the theory of classification of abstract finite-dimensional Lie algebras. 

Historically, Burgers equation is a fundamental partial differential equation from fluid 
mechanics [3]. There have been a number of papers that have contributed to the studies 
of Lie groups of transformations of various classes of Burgers' equation, such as modeling 
of gas dynamics and traffic flow [TT] . 

The classical form of Burgers' equation is 

U t + UU x = VU XX) (1) 

for viscosity coefficient v, which describes the evolution of the field u = u(t, x) under 
nonlinear advection and linear dissipation [3, 10]. When v = 0, Burgers' equation becomes 
the inviscid Burgers' equation: 

u t + uu x = 0. (2) 
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which is a prototype for equations for which the solution can develop discontinuities (shock 
waves). The generalized non-homogeneous inviscid Burgers equation had introduced in 
[T2] has the form 

ut + g(u)u x = f(u). (3) 

In this study, we extend equation ([3]) and consider the following generalized non- 
homogeneous inviscid Burgers' equation (NIB) in the general form 

NIB: u t + g(x,u)u x = f(x,u), (4) 

where / ^ and g ^ are sufficiently smooth nonconstant functions of variables x, u. 

The main results of the present paper are point symmetry and equivalence classifica- 
tion of NIB equations leading to a number of new interesting nonlinear invariant models 
associated to non-trivial invariance algebras. A complete list of these models are given 
for a finite-dimensional subalgebra of infinite-dimensional equivalence algebra derived 
for NIB equations. Reaching to these goals we follow an algorithm explained and per- 
formed in references [H [10]. Our method is similar to the way of [6] for the nonlinear heat 
conductivity equation u t = [E(x,u)u x ) x + H(x,u). 

Throughout this paper we assume that each index of a function implies the derivation 
of the function with respect to it, unless specially stated otherwise. 



2 Principal Lie algebra 

The symmetry approach to the classification of admissible partial differential equations 
relies heavily upon the availability of a constructive way of describing transformation 
groups leaving invariant the form of a given partial differential equation. This is done 
via the well-known infinitesimal method developed by Sophus Lie [8j [9j [10]. Given a 
partial differential equation, the problem of investigating its maximal (in some sense) Lie 
invariance group reduces to solving an over-determined system of linear partial differential 
equations, called the determining equations. 

Considering the total space E = X x U with local coordinate (t, x, u) of independent 
variables (t, x) £ X and dependent variable u £ U, the solution space of Eq. (TjJ, is a 
subvariety S*a C J 1 (K 2 ,M) of the first order jet bundle of 2-dimensional submanifolds of 
E. The one-parameter Lie group of infinitesimal transformations on E is as follows 

t = t + s£(t,x,u) + 0(s 2 ), 

x = x + sr(t,x,u) + 0(s 2 ), (5) 
u = u + s <f(t, x, u) + 0(s 2 ), , 

where s is the group parameter and £, r and tp are the (sufficiently smooth) infinitesimals 
of the transformations for the independent and dependent variables, resp. Thus the 
corresponding infinitesimal generators have the following general form 

d 

v = f (t, x, u) — + r(t, x,u) — + ip(t, x,u)—. (6) 
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According to the invariance condition (see e.g. [TJ, Theorem 2.36), v is a point transfor- 
mation of Eq. if and only if 



u t + gu x -f\ Eq m = 0. (7) 



In the latter relation, is the first order prolongation of vector field v: 

v (1) = V + X, u, u t , u x , u u , u tx , u xx ) + <p x (t, X, u, u t , u x , u tt , u tx , u xx ) — 



I f lu, UJ£, ui X , 'J'ttl <^lXl "'XXI r\ 1 V w / 

OUt ou x 
in which 

^ = B t {Q)+^u tt + ru tx , (9) 
(pF = is x (Q) +£u tx + Tu xx , (10) 

and Q = ip — £u t — r]u x is the characteristic of vector field v and the operators 

= dj + u t ,jdu t + u XyJ du x H 

(for J = t,x) are total derivatives with respect to t and x. Applying v* 1 ) on Eq. @ and 
introducing the relation u t = f — g u x to eliminate u t , we obtain the following determining 
equation: 

-t g g x u x + t f g x - t f x + (p g u u x - (p f u + (p t + f (p u - g (p u u x + g $, t u x - f £ t 
-r t u x - f 2 £ u + f g£ u u x - r u f u x + gip x + gip u u x - gr x u x + g 2 £ x u x - f g£ x = 0. 

In the last relation, since £, r and ip do not depend to variable u x and its powers, so 
Eq. (11 ip is hold if and only if the following equations are fulfilled 

-t g g x + h> g u + 6 g - n + f g£ u - t u f - gr x + g 2 i x = 0, 
Vt-Tf x -(pf u + (p u f-itf-£,uf + gvx-fgix = Q- 

But using the fact that {l,g,g 2 } and {1,/, f 2 } for nonconstant functions f,g are inde- 
pendent sets (see Lemma 1 of [7]), one can divide Eqs. ( IT2|) to the below over-determined 
system 

Tu = 0, £« = 0, Ct-^-T^^O, ipg u -T t = 0, ^ 
£ x = 0, <p x = 0, (f>t - r f x - ip f u = 0, = 0. 

General solution to this system results in the final form of point infinitesimal generators: 

Theorem 1. Infinitesimal generator of every one-parameter Lie group of point symme- 
tries of NIB equation is 

v = |- < 14 > 

Furthermore, every NIB equation in the form Eq. can not be reduced into an inhomo- 
geneous form of a linear equation. 
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Table 1: The commutators table for point symmetry algebra of IBE. 





Vl 


v 2 


v 3 


v 4 


Vl 

v 2 
v 3 
v 4 




V 1 

V F 1 F 2 




V 1 

V F 1 F 2 



V p2 p3 
V 2 

V p2 p4 




-V 3 

V p2 p3 



-V 3 

V p3 p4 


-V 2 
V 3 





Proof: The second assertion is a simple result of p. 209 of [8]. 







Let we change the conditions on /, g and permit some of them be zero: Assume 
that / = and g = g(u). In this case, NIB equation reduced to the the equation 
IBE : Ut + g(u) u x = which has recently studied in [7]. Applying the conditions on /, g 
on determining equation (TTTj) and the same computation as above show that the general 
form of infinitesimal, represented in Theorem 1 of [TJ is not correct. The correct form of 
infinitesimal generators of IBE is as follows 

v = [F 2 (u) t + F\u)}^ + [F\u) g u t + F\u)x + **(«)] + F\u) |-, (15) 

where F % s are arbitrary smooth functions of u. One may divides v to the following vector 
fields 



vi : = 



v 3 := 



'pi 



< F 3 



FHu) 



F A (u) 



d_ 

d_ 

dx ' 



v 2 :-- 



v 4 := 



F 2 



'pi 



F\u) 
- F A (u 



t 9 9 
dt dx 

d d 
ox ou 



(16) 



The commutative table of these vector fields is given in Table 1. 

Also every projective infinitesimal generator introduced in Theorem 2 of [7J has the 
following correct form 



v = [c l( + F'( M )]| + [c 3t + c 1 , + c 2 ]A + ^|., 
for g u 7^ 0. For g(u) = const . , it has the following modified form 



(17) 



v=[cit + F^)]| + [ci, + c 2 ]A + F 2 H |.. 



;is) 



A generalized case of IBE occurs when / = and g = g(x,u). Thus we deal with the 
homogeneous inviscid Burgers' equation. The similar computations show that in this case, 
the general form of point infinitesimal generators is 



v=[F 2 ( M )t + F 1 H]| + e^ 



e 9 dx + F 3 (u) 



d_ 

dx 



(19) 
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3 Equivalence transformations 



It is well-known that an equivalence transformation is a non-degenerate change of the 
variables t, x, u taking any equation of the form (j3J) into an equation of the same form, 
generally speaking, with different f(x,u) and g(x,u): The set of all equivalence trans- 
formations forms an equivalence group E: We shall find a continuous subgroup £ c of it 
making use of the infinitesimal method 0, [TO] . 

To find the group £ c we need to determine those infinitesimal generators 

d d d d d 

Y = f (t, x, u)— + r(t, x,u)— + <p(t, x,u)— + x{t, x, u, f,g)gj + r)(t, x, u, f, g) — , (20) 

from the invariance conditions of Eq. fl4]) as the following system 

J u t + g(x, u) u x = f(x, u) , , 

1 /* = 9t = 0. {2l) 

Here, u and /, g are considered as differential variables: u on the space (t, x) and / , g on 
the extended space (t, x, u). The coordinates ^, r, ip are sought as functions of t, x, u while 
the coordinates x, rj are sought as functions of t, x, u, f, g. The invariance conditions of 
the system (j2T|) are 

/ Y[ut + g(x, u) u x - f(x, u)] = 0, . . 

\ y[/ t ] = Y[g t ] = 0, 

where Y is the prolongation of (|20|) to the first order jet space of differential variables 
t,x,u, f,g which one may represent as follows 

f = Y + V t ^- + V X ^- + X t ^ T + r 1 t ^-. (23) 
du t du x dft og t 

In this relation if 1 and ip x are the same with those introduced in section 2 and 

X 1 = V t ( X ) - ftVm - f x V t {r) - f u V t (cp) = V t ( X ) - f x V t (r) - f u V t (<p), (24) 
rf = V t (rj) - gtV t (£) - g x V u {r) - g u V u (ip) = V t ( V ) - g x V t (r) - g u V u (<p), (25) 

where we consider 

^ d „ d d d j « 

D ' : =a + / 'a/ +£ "95 = a- (26) 

Substituting relations (|24|) - (T26|) in ( j2"3"j) and then applying it on the last two equations of 
(f22l) we find that 

Xt ~ fxT t - f v (p t = 0, (27) 
Vt ~ 9x r t - g u ^t = 0. (28) 

Since the latter equations are hold for every / and g, so we lead to the following facts 

r t = 0, ip t = 0, X t = 0, Vt = 0. (29) 
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Table 2: The commutators table for equivalence symmetry algebra of NIB. 





Y 1 


Y 2 


Y 3 


Y 2 
Y 3 




Y l 


-Y 1 


Y 2 — Y 3 


-Y 1 

Y 3 — Y 2 




By effecting fl23l) on the first equation of (1221 . we have 

-X + VUx + ^ + gtp* = 0. (30) 

Substituting ip t ,Lp x , introducing ut = —gu x + / to eliminate u t and using the fact that 
in the derived relation u x and its powers are free variables, finally we tend to the below 
system 

/ 9 Ca + 9 6 - / r u - g r x + g 2 + rj = 0, (31) 

-x + f<Pu - ftt - f tu + g<f x - fg& = 0. (32) 

The general solution to Eqs. (1231), ([23]) and ([23]) is 

( = £(t,x t u), rj = r(x,u), ip = <p(x,u) (33) 
x = g<px- fgtix- f£t + fPu- f£u, (34) 
v = f T u + gr x - g&- g 2 ix - f g f«. (35) 

After utilizing these relations in Y, one may divide Y to the following vector fields 

r\ r\ r\ 

Y i -= Y i = C-Q t -f[9^ + ^ + fQgj-g[f^u + Ct + gQ-^, 

Y2 --= Y r=^ + {9r x + fr u ]^-, (36) 

Y3 '- =Y ^ = ^o^ + ^^ + f ^df- 

The commutators table of Yi s is given in Table 2. 



4 Preliminary group classification 

In diverse applications of symmetry analysis, the equivalence symmetry is handled to ex- 
tend the principal Lie algebra to the equivalence algebra Lj when some further equations 
under considerations are taken. These extensions are called J-"-extensions of the principal 
Lie algebra. The classification of all non-equivalent equations (with respect to a given 
equivalence group Gj) admitting J-"-extensions of the principal Lie algebra is called a 
preliminary group classification. 

In the general, Gjr is not necessarily the largest equivalence group. In the following, 
we consider a subgroup of the group of all equivalence transformations that has a finite- 
dimensional subalgebra (desirably as large as possible) of an infinite-dimensional algebra 
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Table 3: Commutators table for £10 





X! 


x 2 


X 3 


x 4 


x 5 


x 6 


x 7 


x 8 


x 9 


Xio 













X! 




















x 2 














Xl 





x 2 





x 3 





x 3 

















Xi 





x 2 





x 3 


x 4 


-x x 











-X 4 


-x 5 














x 5 





-x l 





x 4 








-Xi 


-x 5 








X 6 










x 5 














-Xi 


-x 5 


x 7 





-x 2 








x 4 








-x 7 


Xg 





X s 








-x 2 





x 5 





x 7 





Xio — Xq 


x 7 


x 9 





-x 3 











Xi 


— Xg 


Xq — x 10 





x 9 


X\o 








-x 3 








x 5 





-x 7 


— Xg 






with basis f l36j) and use it for a preliminary group classification. We select the subalgebra 
£10, generated by the following vector fields 

d d d 

Xi = — , X 2 = — , X 3 = — , 
at ox ou 

v ^ d , 9 9 v d £ d 2 d 

x * = t di- f df~ 9 W x " = x di' fg df' 9 w 

d d d d d 

Xg = u 7u ~ F 777 ~ f #7T' ^ 7 = x 77" + #7}"' ( 3T ) 
at oj og ox og 

d . d d d 

X 8 = u— + } — , X 9 = 2 — + , 
ax ag aw a j 

Xio = u— + f — , 
aw a/ 

The commutator and adjoint representations of £10 are listed in Tables 3 and 4. 

It is well known that the problem of the construction of the optimal system of solutions 
is equivalent to that of the construction of the optimal system of subalgebras [SI [10] • 
Here, we determine a list (an optimal system) of conjugacy inequivalent subalgebras with 
the property that any other subalgebra is equivalent to a unique member of the list under 
some element of the adjoint representation i.e. i)Ad(g) f) for some g of a considered Lie 
group. The adjoint action is given by the Lie series 

Ad(exp(s Yi))Yj = Yj- 8 [Y h Yj) + j [Y h [Y h Yj]) (38) 



Then we will deal with the construction of the optimal system of subalgebras of C 



10- 
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Table 4: Adjoint table for C 



10 





Xi 


x 2 


x 3 


x 4 x 5 


Xq X 7 


X 8 


Xg 


^10 


A l 


Xi 


x 2 


3 


A A — A 1 A K 


JLq A 7 




Y*_ 
A 9 


A 10 


x 2 


Xi 


x 2 


x 3 


"4 ^5 — sX\ 


X 6 X7 — sX 2 


x 8 


Xq — sX 3 




x 3 


Xi 


x 2 


X 3 


X4 X^ 


X 6 — sXi X 7 


X S 


X 9 


X 10 - sX 3 


x 4 


e'Xi X 2 


*3 


X4 X5 + SX4 


Xq+sXs — — X 4 ^ 7 


X 8 — aX 2 


Xg 


-^10 


Xs 


Xi 


X 2 +sX 


^3 


e _s X 4 X 5 


Xq X 7 + ( e s -l)X 4 


X 8 + sX 5 


X 9 


Xio 


x 6 


Xi 


x 2 


X 3 + sX 1 


X4 — SX5 A'5 


x 6 x 7 


x 8 


s 2 

Xg +SX4 — X5 


X 10 + sX 5 


x 7 


x x 


e s X 2 


X 3 


X4 A' 5 — 3X4 


Xq X 7 


X 8 +sX 7 


e- s Xg 


^10 


x 8 


Xx 


x 2 


X 3 + sX 2 


X4 e~ s X$ 


Xq e~ s X 7 


x 8 


Xg- S (X 10 -X 6 )+( e - s + S - 


-1)X 7 X 1Q + (e- s -lJX 7 


Xg 


Xi 


X 2 +sX 


3 ^3 


X4 X5 


Xq — SX4 X7 + sXg Xq+s(Xq 


s 2 


(Xg— X 4 ) Xg 


XiQ — sXg 


X±Q 


Xi 


x 2 


e s X 3 


X4 Xq 


Xq—sXq X 7 


X 8 +sX 7 


e s Xg 


X±o 



Theorem 2. An optimal system of one-dimensional Lie subalgebras NIB equation in 
the form is provided by those generated by 



1) 


A 1 = 


x u 


11) 


A 11 


= Ti X 6 + X10 


2) 


A 2 = 


x 2 , 


12) 


A 12 


= rj 1 X 8 + X 10 , 


3) 


A 3 = 


x 3 , 


13) 


A 13 


= X 1 + ?7 2 X 8 + X 10 , 


4) 


A A = 


Xi, 


14) 


A u 


= x 2 + Pi x 4 + x 10 , 


5) 


A 5 = 


x 5 , 


15) 


A 15 


= X 2 + 72 X 6 + Xio, 


6) 


A 6 = 


x 6 , 


16) 


A 16 


= aiX 3 + 7 3 X 6 + X 7 , 


7) 


A 7 = 


X3 + X4, 


17) 


A 17 


= a 2 X 3 + 74 X 6 + X 8 , 


8) 


A s = 


X 3 + X 5 , 


18) 


A 18 


= J 5 Xq + T] 3 X 8 + Xg, 


9) 


A 9 = 


Xz + Xq, 


19) 


A 19 


= X4 + 774 X§ + Xio, 


10) 


A 10 = 


- X4 + X10, 


20) 


A 20 


= 7 6 Xq + 775 X 8 + XiQ 



(39) 



where cti, j3i, 7fc, rji (1 < i < 2, 1 < A; < 6, 1 < I < 5) are arbitrary constants. 

Proof: Consider the symmetry algebra £10 of Eq. (J4]), whose adjoint representation was 
determined in Table 4. Given a nonzero vector 

10 

X = Y,a i X i , (40) 

i=i 

we simplify X by eliminating as many of the coefficients a, as possible by use of judicious 
applications of adjoint maps to X. We perform the process through following cases: 

Case a. Let a w 7^ 0. Scaling X if necessary, we can assume that aw = 1: 

X = j2a t X t + X 10 . (41) 

8=1 

When a 3 = or when a 3 7^ by applying Ad(exp(a 3 X 3 ) on X to cancel the coefficient 
of X 3 , one can reduces X to 

X' = aiX 1 + a 2 X 2 + a 4 X4 + --- + Xio. (42) 

Now if 05 = or if as 7^ by effecting Ad(exp(— as X 5 ) on X', we can make the coefficient 
of X 5 vanish: 

X" = ai X 1 + a 2 X 2 + a 4 X 4 + a 6 X 6 + ■ ■ • + X 10 . (43) 
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We can also cancel the coefficient of X 7 . In the case which it is nonzero, for a 7 > 1 by 
applying Ad(exp(ln(a 7 — 1) X 7 ), for a 7 = 1 by applying Ad(exp(— In (2) X 7 ) and for a 7 < 1 
by applying Ad(exp(— ln(l — a 7 ) X 7 ) on X" we change this coefficient to zero: 

X'" = a 1 X 1 + a 2 X 2 + a 4 X 4 + a 6 X 6 + a 8 X 8 + a 9 X 9 + X w . (44) 

Furthermore, for different values of ag, when it is either zero or nonzero, this coefficient 
can be vanished; when ag ^ we act Ad(exp(ag Ag) on X'" to eliminate ag. 

Case al. Let a 8 ^ the additional action of Ad(exp(^ X 4 ) change X to the following 
form which the coefficient of a 2 is deleted 

X"" = a 1 X 1 + a 4 A 4 + a 6 X e + a 8 X 8 + a 9 X 9 + X 10 . (45) 

Case al-1. For a 6 / we effect Ad(exp(^A 3 ). Then the action of Ad(exp(^X 9 ) 
for a 4 ^ or when a 4 = we lead to the following form 

a 6 X 6 + a 8 X 8 + X 10 . (46) 

Now, further simplification is impossible and every one-dimensional subalgebra generated 
by an X with aio, a 8 , a§ ^ is equivalent to the subalgebra spanned by (j47j). this vector 
field was introduced in part 20 of the theorem. 

Case al-2. If we change the condition of Case al-1 to a 6 = 0, the action of Ad(exp(^- Xi) 
for a 4 7^ and then the action of Ad(exp(— hi(^) A 5 ) change X"" to the below final form 
(part 19 of the theorem) 

a 4 X A + a 8 X 8 + X w . (47) 

In this case if a 4 = we can effect Ad(exp(ln(^-) X 4 ) for a\ ^ to reach to part 13 of the 
theorem, while ai = suggests part 12. 
Case a2. Let a 8 in Case al be zero. 

Case a2-l. Moreover, if a 6 ^ by effecting Ad(exp(^ X 3 ) and then Ad(exp(^X 9 ) 
we tend to the form 

a 2 X 2 + a 6 X 6 + X 10 . (48) 

In this form when a 2 7^ the adjoint action of Ad(exp(ln(^-) X 7 ) and when a 2 = we 
find sections 15 and 11 of the theorem resp. 

Case a2-l. Let a$ = 0. Then the coefficient a\ for either a\ ^ by applying 
Ad(exp(ln(^-) X 7 ) or a\ = will be zero: 

X = a 2 X 2 + a A X 4 + X l0 . (49) 

Furthermore, the application of Ad(exp(ln(^-) X 7 ) for a 2 ^ introduces the vector field 
A 2 +a 4 X a +xiq, while for a 2 = by effecting Ad(exp(— hi(^-) X 5 ) to change the coefficient 
of X A equal to 1, we lead to the case A 4 + X w . 

The remaining one-dimensional subalgebras are spanned by vectors of Case a, where 
a 10 = 0. Let Ff(X) = Ad(exp(s X^X) be a linear map, for i — 1, • • • , 10 and every 
X G £10 ■ We continue the classification by a similar method as above. 
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Case b. Let ag 7^ 0, we scale to make ag = 1. One can set the coefficients of 
X 3 , X 2 , X 5 , X 7 , X 1} X4 as zero (when each of the coefficient is zero we eliminate it from 
the list) by effecting F% 3 , F 4 fl2 , F 5 " a5 , Ff a7 , F 3 ai , F 9 a4 resp. Thus we tend to section 18. 

Case c. For a\o = ag = 0,a 8 7^ and assuming a§ = 1, we can make the coef- 
ficients of X 2 , X5, X-?, X\, X4 vanish (when those coefficients are not zero) by applying 
Ft\ F 5 " a5 , F 5 " a5 , Ff a7 , F? 1 , resp. This results in section 17. 

Case d. Let aio = ag = as = 0, a.7 7^ and scale if is necessary to have a 7 = 1, we 
can cancel the coefficients of X 2 , Xl, A 4 , X 5 by actions of F 2 2 , F3 1 , Fg 4 , F^ resp. to have 
section 16. 

Case e. Let aio = ■ • • = a 7 = 0, a 6 7^ 0. We can suppose that a 6 = 1. Applying 
F 3 ai , Fg 4 , Ffg , Fg a2 on X we can changed the coefficients of X 1; X 4 , X 5 , X 2 to zero (if they 
are not zero). Then if 03 7^ by acting Ad(exp(ln(^) X w ) we find section 9 and when 
a 3 = we lead to section 6. 

Case f. Suppose that aio = • • • = ^6 = 0, as 7^ 0. By assuming as = 1, we can 
make the coefficients of X\, X4, X 2 vanish by F 2 l , -F 4 _a4 , Fg a2 . Then if 03 7^ we apply 
Ad(exp(ln(^-) X 6 ) to give section 8 and for a 3 = we have section 5. 

Case g. Consider the case which aio — • " — a 5 = 0, a 4 7^ 0. We can assume that 
a 4 = 1. we can make the coefficients of Xx, X 2 vanish by F^ 1 , F% a2 . In addition for a 3 7^ 0, 
scaling the coefficient of A 3 by Ad(exp(ln(^) X 10 ) we have section 7, while 0,3 = leads 
to section 4. 

Case h. For a 10 = ■ ■ • = a 4 = 0, a 3 7^ and assuming a 3 = 1, the coefficients of Xi, X 2 
will be vanished by the actions of F 6 ~ ai , F 8 ~ a2 and X reduces to section 3. 

Case i. If aio — • • • — a 3 — 0, a 2 7^ 0. Scaling a 2 = 1, by effecting F§ ai the coefficient 
of X\ will be canceled which suggests section 2. 

Case j. Finally when aio = • • • = a 2 = we have section 1 of the theorem. 

There is not any more possible case for studying and the proof is complete. <) 

Now, since /, g in Eq. (Tj0) are functions of variables x, u, so we project their derived 
optimal system to the space (x,u,f,g). The nonzero vector fields of (x, u)— projections 
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of dSSJ follows 



1) 


Z 1 = 


= A 2 = 


d x , 


2) 


Z 2 -- 


= A 3 = 


ud u , 


3) 


Z 3 -- 


= A 4 = 


fdf + gdg, 


4) 


Z 4 -- 


= A 5 = 


g[fd f + gd g ], 


5) 


Z 5 -- 


= A 6 = 


f[fd f + gd g ], 


6) 


Z 6 -- 


= A 7 = 


d u -fd f -gd g , 


7) 


z 7 = 


= A 8 = 


d u - g[f d f + g d g ], 


8) 


z 8 = 


= A 9 = 


du-fifd f + gdg}, 


9) 


z 9 = 


= A 10 = 


= ud u - gd g , 


10) 


z 10 


= A u 


= ud u + (l -7i f)fd f - f gdg, 


11) 


z 11 


= A 12 


= A 13 = r] 1 ud x + ud u + fdf + ^gdg 


12) 


z 12 


= A 14 


= d x +ud u + (1 - Pt)f d f - fiigdg, 


13) 


z 13 


= A 15 


= d x + u d u + /[(l - j2f)d f -j 2 g d g ] , 
= xd x --f 3 f 2 d f + (l- 7s f)gd g , 
= ud x + a 2 d u - 74 f 2 d f + (1 - 74 g)f d g , 
= rj 3 u d x + (x + 75 u) d u + (g- 75 P)d g , 


14) 


z 14 


= A 16 


15) 


z 15 


= A 17 


16) 


z 16 


= A 1S 


17) 


z 17 


= A 19 


= r] i ud x + ud u + (774 / - flf)5 fl , 


18) 


z 18 


= A 20 


= Tj^u d x -\- u d u + /[(l - 7 6 /) 9/ + (775 - 



(50) 



According to paper 7 of [3j we conclude that 



Proposition 3. Let C m :— (Xi : % — 1, • • • ,m) be an m- dimensional algebra. Denote 
by A 1 (i = 1, • ■ ■ , s, < s < m, s 6 N) an optimal system of one-dimensional subalgebras 
of C m and by Z % (i = 1, • • • ,t, < t < s, t G N) t/ie projections of A 1 , i.e., Z % = pi(A l ). 
If equations 

g = g(x,u), f = f(x,u), (51) 
are invariant with respect to the optimal system Z % then the equation 

u t + g(x,u) u x = f(x,u), (52) 
admits the operators Y % = projection of A 1 on (t,x,u). 



Proposition 4. Let Eq. / f5^j) and the equation 

u t + g(x,u)u x = f(x,u), (53) 

be constructed according to Proposition 3 via optimal systems Z % and Z resp. If the subal- 
gebras spanned on the optimal systems Z l and Z* resp. are similar in C m , then Eqs. ( T5T3I) 
and (53}) are equivalent with respect to the equivalence group G m generated by C m . 



Regarding to Propositions 3 and 4, now our task is the investigation for all non- 
equivalent equations in the form of Eq. (j3J) which admit ^-extensions of the principal Lie 
algebra Cg, by one dimension, that are equations of the form (TJJ such that they admit, 
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together with the one basic operator (|14p of C±, also a second operator X^ 2 \ In each 
case which this extension occurs, we indicate the corresponding coefficients /, g and the 
additional operator X^ 2 \ 

By using the algorithm for operators Z % {% = 1, • • • , 18) to /, g one can find vector fields 
s. Let we consider the following examples. 

We select the operator 



Z 5 = f 2 4j + f9-^- (54) 
df dg 



The characteristic equations are 



P ~ fg' { j 



so the invariants of Z are 



/ 

h = x, I 2 = u, h = -■ (56) 

9 

In this case there are no invariant equations because the necessary condition for existence 
of invariant solutions is not fulfilled (see [TO], Section 19.3), that is, invariants (156 \ cannot 
be solved with respect to / and g since I 2 is not an invariant function of I\. Similar results 
are hold for vector fields Z 3 and Z 4 . 

Let perform the algorithm for another example, by considering vector filed 

zu =4,-^ 2 w + ^K' (57) 

then the characteristic equations corresponding to Z 14 is 

— = d f = dg (w\ 

which determines invariants for 73 7^ 0. Invariants can be taken in the following form 

h = u, / 3 = /-_?_, h = ^-j. (59) 
73 mx x l 73 1 

From the invariance equations we can write 

/ a = $(/i), J 3 = *(A). (60) 

It results in the forms 

f = <&(\) + —L- g = x 1 ~ 13[HX)+ ^ ] m(\), (61) 

73 mx 

where X = u. For the case which 73 = 0, Z14 is in the form xd x + gd g and similar way 
shows that / = $(A) and g = x\&(A) when X = u. 
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Table 5: The result of the classification 



N Z 



Invariant A Equation 



Additional operator 



1 
2 
3 
4 
5 
6 
7 



z 2 
z 6 
z 7 
z 8 
z 9 

z io 

z 11 



9 Z 12 

10 Z 13 ( 72 

11 Z 13 ( 72 

12 Z 14 ( 73 

13 Z 14 ( 7 = 

14 Z 15 (q 2 

15 Z 15 (q 2 

16 Z 15 (a 2 

17 Z 15 (q 2 

18 Z 16 (t? 3 

19 Z 16 ( m 

20 Z 16 (t? 3 

21 Z 16 (t? 3 

22 Z 17 (t? 4 

23 Z 17 (t? 4 

24 Z 18 ( % 

25 Z 18 ( % 

26 Z 18 ( % 

27 Z 18 ( % 



7^0) 
= 0) 

7^0) 
= 0) 

= 74 = 0) 

7^0,74 = 0) 

= 0, 74 ^0) 

/0, 7 4y^0) 

= 75 = 0) 
7^0,75 = 0) 

= 0,75 7^0) 
7^0,75^0) 
7^0) 
= 0) 

= 76 = 0) 
7^0,76 = 0) 
= 0,76 7^0) 
7^ 0,76 7^0) 



u 

X 
X 
X 
X 
X 
X 

x — rj±u 

Inu — x 

Inu — x 
lnu — x 



u t + ^u x = <& 
u t + ^>u x = <$> 



u t +[#+l/ii]u a . = [¥+l/u]$ 
u t + [$ + 1/u] $>u x = $ + l/u 



u t +u^u x = & 



JL 

dx 

d_ 

du 

1 dt ^ Qu 

r IL , d_ 

x dt ~ du 

" dt ^ du 

+ JL -l 7 , JL 

dt 



9m 



^±(7i/-l)*^- w ^ 
lit + -r]iu<&\u x = u$ 

u t + ti /3l ^'ti :r = n 1 ~' 3l $ 
^+[(^1-1)*]-^ = 



72 
72 u$— 1 



73 In a; 



72«<E>— 1 

^+[7^(^-1)*]^=*-^ 

u t +^[l±$*]n a; = $ 
ut + [^+e u ' x ]u x = <S> 
±(r] 3 u 2 -x 2 ) u t + l^+e u / x ]u x = § 

u t + [~ 1b <$> 2 ±{~i b u+x)~^]u x = § 
u t +[ 75 $ 2 ±e a; /( ? ? 3U )+*] u x = $ 
u t + [rj4_<&±u^>\u x = <& 



u 
a 2 



U 



x is 



JL 

dx 

d_ 

dx 



+ a 2 



On 



a 2 



2 
X 

B 



X 



- u 



Inu ± x 

x 



u— - 



V5 



74^| + u £ 

x^- 

x du 

^ u m + x ii 

hsu + x}-^ 

V3U-^ + [7 5 U + x}-^ 

L dt ^ a du 
U du 

^ U m + u ii 
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From Proposition 3 applied to the operator Z 14 we obtain the following additional 
operators for the introduced functions / and g as above 

* <2, ^4+4- < 62 > 

When 73 7^ 0, while for 73 = we find = xd x as an additional operator. 

Repeating the algorithm for other Z % s of fl50l) . the preliminary group classification of 
an NIB equation of the form (j3J) admitting an extension £2 of the principal Lie algebra 
L\ is listed in Table 5. In this table we assume that 

B = x + j 3 u + — LambertW(A) (63) 

73 



where 

^ 3 e -Tf/»n 

T]iU 

and $ = $(A), ^ = ^(A) are arbitrary functions of the invariant A 



A = — iLe-ni/n, (64) 



5 Conclusion 

A symmetry analysis for generalized non-homogeneous inviscid Burger's equations of 
class (j3J) led to find the structure of point infinitesimal generators as well as equivalence 
operators. In addition we found the modified forms of point and projective symmetries of 
IBE in the paper [7] among with the general form of point infinitesimals of a generalized 
case with respect to IBE and results of p2]. Also, equivalence classification is given of 
the equation admitting an extension by one of the principal Lie algebra of the equation. 
The paper is one of few applications of a new algebraic approach to the problem of 
group classification: the method of preliminary group classification. Derived results are 
summarized in Table 5. 
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